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The problem of the effects of compressibility and large-scale anisotropy on anomalous scaling 
behavior is considered for two models describing passive advection of scalar density and tracer 
fields. The advecting velocity field is Gaussian, (5-correlated in time, and scales with a positive 
exponent e. Explicit inertial-range expressions for the scalar correlation functions are obtained; 
they are represented by superpositions of power laws with nonuniversal amplitudes and universal 
(dependent only on e and a, the compressibility parameter) anomalous exponents. The complete set 
of anomalous exponents for the pair correlation functions is found nonperturbatively, in any space 
^ , dimension d, using the zero-mode technique. For higher-order correlation functions, the anomalous 

' exponents are calculated to O(e^) using the renormalization group. Like in the incompressible case, 

the exponents exhibit a hierarchy related to the degree of anisotropy: the leading contributions to 
the even correlation functions are given by the exponents from the isotropic shell, in agreement with 
the idea of restored small-scale isotropy. As the degree of compressibility increases, the corrections 
become closer to the leading terms. The small-scale anisotropy reveals itself in the odd ratios of 
correlation functions: the skewness factor is slowly decreasing going down to small scales for the 
incompressible case, but becomes increasing if a is large enough. The higher odd dimensionless 
ratios (hyperskewness etc.) increase, thus signalling the persistent small-scale anisotropy; this effect 
. becomes more pronounced for larger values of a. 

PACS number(s): 47.27.Te, 47.27.-i, 05.10.Cc 
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^ . I. INTRODUCTION 

<N 

' Much attention has been paid recently to a simple model of the passive scalar advection by a self-similar Gaussian 
' white-in-time velocity field, the so-called "rapid-change model," introduced by Kraichnan see, e.g., Refs. (2|-||] and 
' references therein. Despite its simplicity, the model reproduces many of the anomalous features of genuine turbulent 
' heat or mass transport observed in experiments. On the other hand, it appears easier tractable theoretically: for 
<^ . the first time, the anomalous exponents have been calculated on the basis of a microscopic model and within regular 

^ ■ expansions in formal small parameters Therefore, the passive advection by the "synthetic" velocity with 

. prescribed statistics, being of practical importance in itself, may also be viewed as the starting point in studying 

C ■ anomalous scaling in the turbulence on the whole. 

^ , In the original Kraichnan model, the velocity field is taken to be Gaussian, isotropic, incompressible and decorrelated 
'k^ in time. More realistic models should involve anisotropy and compressibility. Recent studies have pointed out 
rN significant differences between the compressible and incompressible cases |l0|-|l7| . It is noteworthy that the potential 
■ velocity field remains nontrivial in the one-dimensional case, which is more accessible to numerical simulations and 
" " ' allows interesting comparison between the numerical and analytical results; see Ref. |]lo| . 

Another important question recently addressed is the effects of large-scale aniso tropy on inertial-range statistics of 
passively advected scalar p^[-p3| and vector [p4[-p6| fields and the velocity itself |2^j2^. These studies have shown 
that the anisotropy present at large scales has a strong influence on the small-scale statistical properties of the scalar, 
in disagreement with what was expected on the basis of the cascade ideas |l8|-p0|. On the other hand, the exponents 
describing the inertial-range scaling exhibit universality and hierarchy related to the degree of anisotropy, which gives 
some quantitative support to Kolmogorov's hypothesis on the restored local isotropy of the inertial-range turbulence 
0,|1]-|28|. 



In this paper we analyze the effects of the large-scale anisotropy induced by a random source on a passive scalar 
by two methods: First, we carry out the zero-mode calculation of the correlation function of the passive scalar with 
an anisotropic source field and in an isotropic compressible velocity field decorrelated in time. Second, for the same 
model, we have performed the two-loop renormalization-group analysis of the asymptotic behavior of the structure 
functions of the passive scalar of arbitrary order. This paper is organized as follows: In Sec. II, the zero-mode solution 
for the correlation function both for passive density and passive tracer is constructed. Two-loop renormalization-group 
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analysis of the structure functions of the passive scalar is carried out in Sec. Ill with the use of the operator-product 
expansion. Section IV is devoted to the discussion of the results. 



II. ZERO-MODE SOLUTION FOR PASSIVE DENSITY AND TRACER 

There are two types of diffusion-advection problems for the compressible velocity field |^^. Passive advection of a 
density field 0{x) = 0{t, x) (say, the density of an impurity) is described by the equation 

dt9 + d,{v,9)^i^„d^e + f, (1) 

while the advection of a "tracer" (say, temperature, specific entropy, or concentration of the impurity particles) is 
described by 

dte + {uA)d = lyod^e + f. (2) 

Here dt = d/dt, di = d/dxt, vq is the molecular diffusivity coefhcient, is the Laplace operator, v(x) is the velocity 
field, and / = f{x) is an artificial Gaussian scalar noise with zero mean and the covariance 

{f{x).f{x'))=5{t-t')C{v), r^x-x', (3) 

where C(r) varies noticeably on r = |r| ~ L, the integral turbulence scale. In the presence of a preferred direction 
specified by a unit vector n, the function C(r) can be written in the form 



C{v)^Y.^i{mr)Pi{z), z^^, (4) 



where m = \/L, Ci{mr) are coefficient functions such that C(r) becomes constant at mr = and decays rapidly for 
mr ^ cx), z is the cosine of the angle between n and r, and Pi{z) are (d-dimensional) Legendre polynomials satisfying 
the equations 

(1 - z2)P/'(z) + z (1 - d)P/(z) + l{l + d- l)Piiz) = 0. (5) 

The anisotropy makes it possible to introduce also a mixed correlator (v/) cx n6{t — t') C"(r) with some function 
C similar to C in Eq. (^. This violates the evenness in n and gives rise to nonvanishing odd correlation functions 
of 9, but leads to no serious alterations in the analysis. We shall discuss this case later on, and for the time being we 
assume (v/) = 0. 

In the real problem, the field v(x) satisfies the Navier-Stokes equation. In the simplified model considered in 
it obeys a Gaussian distribution with zero mean and the covariance 



with 



{v,{x)v,{x')) ^5{t-t')K,j{r) (6) 

where /ij(k) = Sij — kikj/k'^ and Qij(k.) — k^kj/k'^ are the transverse and longitudinal projectors, respectively, 
k = |k|. Do and Dq are positive amplitude factors, and d is the dimensionality of the coordinate space. For D'q = 
(the incompressible case) the models (|l|) and (^ coincide. For < e < 2, the so-called eddy diffusivity 



S,,{r)^K,,{0)~K,,iv) (8) 



has a finite limit for 



with 



S,,{r)=Dr' {d + e-l + a)d,j-e{a-l)^ , (9) 



^=(4.)^/22^(d + e)r(d/2 + s/2)' " = ^o/^°' (1°) 
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where T is the Euler gamma function (note that D and a both are positive). In the renormahzation group (RG) 
approach, the exponent e plays the same role as the parameter e = A — d does in the RG theory of critical behavior; 
see Refs. 16|. The relation Dq/vq = defines the characteristic ultraviolet wave-number scale A. 

The issue of interest is the behavior of various correlation functions in the incrtial range specified by the inequalities 
Ar ^ 1, mr ^ 1. In the models odd multipoint correlation functions of the scalar field vanish, while the even 

equal-time functions satisfy linear partial differential equations; see, e.g., The equation for the equal-time pair 

correlation function D{r) = {9{t,x)9(t,x')) , is easily derived from the Schwinger-Dyson equations (see Refs. jr6| , p5| ) 
and has the form (here and below in equal-time functions, we omit common to all the quantities time arguments): 

2j,od^D{r) + [S,,{r)d,dj] D{v) = C(r) (11) 

for the model (||) and 

2i,od^D{r) + d^dj [S^j{r)D{v)] = C(r) (12) 

for the model (|l|), with C(r) from (^ and Sij{r) from (|^). In the presence of the preferred direction, n, the correlation 
function can be decomposed in the Legendrc polynomials, 

oo 

Z?(r) = ^A(r)P,(^), (13) 
/=o 

where the coefficient functions are sought in the powerlike form 

Di{r)c^Dir^K (14) 

Owing to the evenness in n, only even polynomials contribute to ([l^). 

It is well-known (see, e.g., Refs. [pH6l) that the nontrivial inertial-range exponents are determined by the zero 
modes, i.e., solutions of Eqs. (|ll|), (^2^ neglecting both the forcing [C(r) = 0] and the dissipation (i/q = 0). The 
homogeneous equations are SO{d) covariant, and the equations for the coefficient functions in ([l^) foliate. Substituting 
the representations (p^), (|lj) into Eqs. (p]), (|l^ ) and using the relations diSij{r) — ae{d + £)Drjr^'^^^ , didjSij{r) = 
ae{d + e){d — 2 + £)Dr^'^^^ then gives quadratic equations for the exponents C,i in Eq. (p^), namely, 

0(0+d-2)-U^ + rf-2) + ^^^f-iM^-^=0 (15) 

[a — 1 + a + e) 

for the tracer, which has two solutions, — 2 — d — I ^ 0{e) and 

lil - l){a - 1) 9, 
^'-^-^ (d + L2)(d-l + .) +^^^)' ^''^ 

and 

0(0+c^-2)-/(^ + .-2)+e ^'(^'-^)(-;) + "^-^ + ^^^^'+-^-^ + -) ^0 (17) 

[a — \ + a + £) 

for the density, with two solutions, Qi = —d — I + 0{e) and 

Jil-l){a-l) + adid + 2l-2) 
id + 2l-2)id-l + a) + (^^^ 

The standard arguments f^-j^ show that only the second solution, 0=1 + 0{e), is "admissible." It has the form 

C; = [2(-l + d + a + ae)]"^ {-2 + 3d - d'^ + 2a - da + e - de + ae 

+ [{2-3d + d'^ -2a + da-e + de-aef (19) 
~Ali~2 + 3d-d'^ + l-dl + 2a-da-la + 2e -de-l£){-l + d + a + ae)]^^'^^ 

for the tracer and 
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[2 (-1 + d + a + ae)]~^ {-2 + 3d - (f + 2a - da + e - de + ae - 2dae - 2ae'^ 
(2 - 3d + d"^ ~ 2a + da - e + de - ae + 2dae + 2ae^y - 4 (-1 + d + a + ae) 

X {-2l + 3dl - d"^ I + 1^ - dl'^ + 2la- dla~ a + 2le - die - f e -2dae (20) 

+d'^ae-2as^ + 2 da + a e^)]^/^^ 

for the density. For a = 0, these sohitions coincide with each other and with the exponents obtained earher in 
Refs. for the incompressible case. 

The exponents (|l^), ( pO| ) exhibit a hierarchy related to the degree of anisotropy: 

CiXi' if l>l', (21) 

i.e., the less is the index I, the less is the exponent and, consequently, the more important is the contribution to the 
inertial-range behavior. The leading term is given by the exponent from the "isotropic shell." This behavior is 
illustrated by Figs. 1 and 2. 

Although the hierarchy holds for all values of a {dQ /dl > 0) , the corrections become closer to leading terms as a 
increases: d^C,i/dlda < 0. This behavior is illustrated by Fig. 3. 

Since the equation (|^) is invariant with respect to the shift 6^9 + const, the relevant quantities for the tracer are 
the so-called structure fmictions, 

5„(r) = (Ki,x)-0(t,x')]"), r = x-x', (22) 

with the Legendre decomposition 

oo 

5„(r) =^5„,r';-Fi(^) (23) 

(=0 

with some numerical coefficients Sni- The comparison with Eqs. (|l3|), ( [l^ ) gives C21 — Ci with Q from ( pj| ) for all 
I > 0. For I = 0, the constant term with Co = drops out from the difference in (p2|), and the behavior of the isotropic 
shell is determined by the sublcading exponent (2 — e); see, e.g., Rcf. [0. Note that the hierarchy relations ( |2l| ) 
remain valid also for 5*2 = [D{0) ~ D{r)]/2. 
For the density case, the exponent 

^ -e{d + e)a 
^" {d-l) + a{l+e) 

(the square root in Eq. ( ^o|) is taken explicitly) gives the leading contribution both for the pair correlation function 
D{r) and the structure function S2 in ([2^), in agreement with the exact solution of ||l6| (for d — I, see Ref. [p0[). 
Note that for this case, the anomalous scaling emerges already for the pair correlation function, like in the model of 
passively advected magnetic field studied in Ref. p9 . 



III. TWO-LOOP RENORMALIZATION-GROUP ANALYSIS OF STRUCTURE FUNCTIONS 

The higher-order structure functions can be studied using the field theoretic renormalization group and oper- 
ator product expan sion. The detailed exposition of these techniques and practical calculations can be found in 
Refs. [p| Jl^Jl7| , |2l[p5|| ; below we confine ourselves to only the necessary information. 

The field theoretic models corresponding to the stochastic equations (^ and (|^) are multiplicatively renormaliz- 
able; the corresponding RG equations have infrared stable fixed points. In particular, this leads to the following 
representations for the structure functions in the model (|^) in the inertial range (Ar ^ 1, mr <C 1): 

^„ = ^a(r,z)(mr)^". (24) 

a 

Here Ca(r, z) are coefficients analytical in m and finite for to — > 0, and are the critical dimensions of the composite 
operators entering in the operator product expansion. 

The leading zero-mode contribution in the Z-th shell for 5„ is determined by the critical dimension A„; of the 
irreducible traceless l-th rank tensor operator built of n fields 9 and minimal possible number of derivatives. For I < n 
such an operator has the form 
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l + 2p. 



(25) 



Here the dots stand for the appropriate subtractions involving the Kronecker delta symbols, which ensure that the 
resulting expressions are traceless with respect to contraction of any given pair of indices, for example, diOdjO — 
Sijdkddkd/d, diOdjOdk.O — {SijdkO + SikdjO + Sjkdi6)/{d + 2) and so on. 

The exponents Cni — "■(! ^ e/2) + A„i are calculated in the form of scries in e, where (ni = n + J2T=i Cni' 
the first order in s: 



1 



(d + 2) 



(n ~ l){d + n + l) 1(1 - l)(a - 1) + a{n - l){n + l- 2) 
2 ^ (d - 1 + a) 



For the incompressible model, the exponent from the isotropic shell [l — 0) was obtained in Refs. 



0{l/d) and in Refs. 
in Ref. |^l[|; see also 
of those papers, p = 



(26) 



to the order 



to the order 0(e); the results for n — 2> are given in Ref. M. The general case can be found 
2^3l[]. For general a > 0, the exponent Cno found in Ref . 1 1^ (see also |l^]; in the notation 
a/ld — 1 + a)). The result for general n, I is given in Ref. [Ell (for more details, see |17[|). 



We have performed the two-loop calculation of the exponents C,ni and obtained: 



a(2) 



'd{d + 1) + (d^ -2d- 4)a + {3d + 4)^2] (d^i - kz) 



d{d + 2)^{d- 
in-2) 



8d[-(d2 + 5d4 
24) - {d + A){d^ 



4d(d + 2)3((i + 4)(d - l + ( 
A[{d + 1){M^ + nd^ + 2M 

3{d + 2)/i(d)|4d[3 + {d~ 4)a - {d - l)a^]Ki + [- 
2{d^ + lOd^ - d - 4)a - 3(2^2 _ 3d - 4)a2]K2|) 



10) - {d-2){d + A)a 

h 7d^ -2d- 4)a + {d 
[-i{d+l){d^ + bd- 



- (2d2 + 
l){bd^ 



ld + 2)a^]ni 
f 8d-24)a2]K2 



(27) 



4) 



where we have written ki = n{n — 1), K2 = (n — l){n + ? — 2 + d), h{d) — F(l, 1, d/2 + 1; 1/4) and F(a, 6; c; z) 
is the hypergeometric function; see, e.g., Ref. |^^. For integer space dimension d one has h{\) — 27r/(3-\/3) and 
h{2) = 41n(4/3), and the others can be obtained from the recursion relation 



ih{d) 



d + 2 



h{d + 2) ^ 4, 



(28) 



valid for all d. 

For the incompressible case, a = 0, Eq. (|2 



.(2 



) _ (d+ l)(dKi - K2) 

~ (d + 2)3(d-l)2 " 

4[(d+l)(3d3 + 17d2 



becomes: 

{n-2) 



Ad{d 
-20d- 



h2)3(d + 4)(d- 
24)]K2 + 9(dH 



-^{-8d(d2 



L) 

2)h{d)[AdKi 



f 5d 

(d4 



f 10)ki 
l)(d2 + 5d- 



4)k2]} 



(29) 



The results for I = and 2 were obtained earlier in Ref. ||^. In that paper, they were expressed in terms of the 
functions hi{d) = h{d + 2) and h2{d) = h{d + 4), which can be reduced to the form (29) using the relation (p^). 

For d = 1, only the exponents C„o and C„i make sense (traceless operators of the rank I > 2 vanish identically). 
They are independent of a (the velocity is purely potential) and have the form (/ = if n is even and Z = 1 if n is 
odd) 



"2 



n{n - \){n - 2)e^n 
6\/3 



(30) 



The expression ( ^0| ) is in agreement with the 0{e) result of Ref. [|oj and the ©(e^) result of Ref. ||T^ (in the latter, 
the density case was studied, but in one dimension these models can be related by the replacement 9 89). 

For I > n, the leading operators contain more derivatives than fields, the corresponding exponents behave as 
Cni = I — n + 0(e) and the corresponding terms in representation (^3|) rapidly decrease for mr 0. 

Now let us return to the density model (|^). It is not invariant with respect to the shift 9 6* + const, the operators 
0" have nontrivial critical dimensions. 



A„ = n (-1 + I) 



an{n — l)d£ a{a — l)n{n - l){d — l)e'^ a'^n{n — l){n — 2)dh{d)e'^ 



2(d- 1 



2(d- 1 



4(d- 1 



(31) 
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and different terms in the structure functions have different scahngs, see Ref. [Q. The relevant quantities are then 
the equal-time pair-correlation functions of the powers of 6, which have the form 

(6l"(x)6'?'(x')) oc v-^'^+^^''^ A-("+P) (Ar)-'^"-'^" Kp(r)(Mr)'^"+''. (32) 
The leading terms in the Legendre decompositions for the scaling functions, 

C30 

i^„,(r) = ^c("''\mr)F,(^), (33) 
;=o 

are given by the contribution from the scalar operators 0"+^ without derivatives, Cq"^-* cx (mr)^"+''. 

The operators that determine corrections (/ > 0) to the leading term with / — necessarily contain I derivatives 
and, in contrast to the tracer case, the correction exponents differ from the leading ones by a few unities, I -f 0(e). 
In particular, the leading / = 2 correction is related to the operator 



whose dimension equals to 

_ _£_ I (liiiiMf^ , fci) } , oi,r 

Note also that the exponents Q in Eqs. (|l^), (^0|) are related to the composite operators with two fields 9 and I 
free indices, which (up to total derivatives and subtractions with the delta symbols) reduce to the form 9di-^ ■ ■ ■ di^9; 
the one-loop calculation confirms the 0{e) results ([l6|), (|l8|). 

Since the leading terms of the even functions ( p3| ) are determined by the exponents of the isotropic shell (i.e., those 
related to scalar composite operators), the inertial-range behavior of the former is the same as in the isotropic model. 
This gives quantitative support both to Kolmogorov's hypothesis on the restored local isotropy of the inertial-range 
turbulence and to the universality of anomalous exponents with respect to the way in which the turbulence is excited. 

On the contrary, the small-scale anisotropy reveals itself in odd correlation functions (which are nonzero in the 
presence of a mixed correlator (v/) or the constant mean gradient of the scalar field). It follows from the above 

analysis that the dimcnsionless ratios TZn = S2n+i/ S'^"'^^^^'^ for the tracer case in the inertial range have the form 

7^„ cx (mr)'^^"+i'S (34) 



where A2„+i,i = C2ri+i,i — (2n+ 1)(2 — e) is the critical dimension of the vector composite operator ( |25| ) built of 2n-|- 1 
scalar gradients. From Eq. ( [26| ) we find, to the first order in e, 

A2„+i,i = e[{d~l + a){d + 2- Aii^) - 8a n^] /2{d + 2){d-l + a). (35) 

For small a, the skewness factor TZi decreases for m,r 0, but slower than expected on the basis of cascade ideas (the 
latter suggest that the odd ratios should vanish for small mr, where the turbulence is expected to become isotropic), 
while the higher-order ratios increase, thus signalling the persistence of small-scale anisotropy. When a increases, TZi 
also becomes divergent for mr provided a is large enough [namely, a > {d — l){d + 2) / {10 — d) + 0(e)], while the 
higher-order ratios diverge even faster. 

It was argued in Refs. ||l3-n5[ that the anomalous scaling regime in the models at hand breaks down if e and a are 
both large enough [p = a/Jd — 1 + a) > d/^'^] and the inverse energy cascade with no anomalous scaling takes place. 
This effect obviously cannot be detected within the e expansion. It is noteworthy that the exact nonperturbative 
exponents (^9|), ( po|) show no hint of anomaly at the threshold, p = (i/^^, in contrast to the exact exponents for the 
magnetic case which become complex when the anomalous scaling regime breaks down; see Refs. pO|pl]. 



IV. CONCLUSION 



To conclude with, we have studied the effects of compressibility and large-scale anisotropy on the anomalous scaling 
behavior in two models, which describe passive advection of scalar tracer and density fields. The advecting velocity 
field is Gaussian, (5-correlated in time, and its spatial correlations scale with a positive exponent e. Explicit inertial- 
range expressions for the scalar correlation functions have been obtained; they are represented by superpositions of 
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power laws with nonuniversal amplitudes and universal (dependent only on e and a, the compressibility parameter) 
anomalous exponents. The complete set of anomalous exponents for the pair correlation functions has been found 
nonperturbatively, in any space dimension d, using the zero-mode approach. For higher-order correlation functions, 
the anomalous exponents have been calculated to 0{e^) using the RG techniques. Like in the incompressible case, the 
exponents exhibit a hierarchy related to the degree of anisotropy; the leading contributions to the even correlation 
functions are given by the exponents from the isotropic shell, in agreement with the idea of restored small-scale 
isotropy. 

This picture seems rather general, being compatible with that established recently for the Navier-Stokes turbulence 
p8[ , passive scalar advection by the two-dimensional Navier-Stokes field j2^, and passive advection of a scalar |2^] 
and vector |^,^ fields by the white-in-time incompressible synthetic velocity field. 

As the degree of compressibility increases, the corrections become closer to the leading terms; cf. Ref. for the 
passive advection of a magnetic field. 

On the contrary, the small-scale anisotropy reveals itself in the odd ratios of correlation functions: the skewness 
factor is slowly decreasing down to small scales for the incompressible case 0, but becomes increasing if a is large 
enough. The higher-order odd dimcnsionless ratios (hyperskewness etc.) increase, thus signalling the persistent small- 

|. This effect becomes even more pronounced for larger values of a; cf. Ref. p^ ] 
for the magnetic case. 



scale anisotropy, cf. Refs. |17^;E1 23 
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FIG. 1. Behavior of the exponents Q (Z = 2, 4 and 6 from the bottom to the top) from Eq. (^^ vs e in three dimensions for 
= (thin lines) and a = co (thick lines). 
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FIG. 2. Behavior of the exponents C,i {I = 0, 2, 4 and 6 from the bottom to the top) from Eq. ( po| ) vs e in three dimensions 
for a = Q (thin lines) and a = oo (thick lines). 
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FIG. 3. Behavior of the exponents ^; {I — 2, A and 6 from the bottom to the top) from Eq. nLm vs p = a/(d — 1 + a) for 
= 1, d = 3 (thin hnes) and d = 2 (thick hnes). 
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